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We present a microscopic laser model for many atoms coupled to a single cavity mode, including the light 
forces resulting from atom-field momentum exchange. Within a semiclassical description, we solve the equa- 
tions for atomic motion and internal dynamics to obtain analytic expressions for the optical potential and friction 
force seen by each atom. When optical gain is maximum at frequencies where the light field extracts kinetic 
energy from the atomic motion, the dynamics combines optical lasing and motional cooling. From the cor- 
responding momentum diffusion coefficient we predict sub-Doppler temperatures in the stationary state. This 
generalizes the theory of cavity enhanced laser cooling to active cavity systems. We identify the gain induced 
reduction of the effective resonator linewidth as key origin for the faster cooling and lower temperatures, which 
implys that a bad cavity with a gain medium can replace a high-Q cavity. In addition, this shows the impor- 
tance of light forces for gas lasers in the low-temperature limit, where atoms can arrange in a periodic pattern 
maximizing gain and counteracting spatial hole burning. Ultimately, in the low temperature limit, such a setup 
should allow to combine optical lasing and atom lasing in single device. 



I. INTRODUCTION 

Light forces in optical resonators are strongly modified as 
compared to free space | ll0|3|. Extending the ideas of laser 
cooling in free space, various new schemes of laser cooling 
of atoms have now been developed in the field of cavity QED 
theoretically (for a review, see e.g. Ref. |4]) and experimen- 
tally |5|. The common basis of all these considerations consti- 
tute the strong correlations between the atomic motion and the 
dynamics of the cavity field. Contrary to conventional laser 
cooling schemes relying on spontaneous emission, dissipa- 
tion of energy and entropy from the atomic motion is provided 
here by the cavity decay channel, where the emission of blue 
shifted photons is enhanced. Under ideal circumstances, it 
is possible to implement very efficient laser cooling schemes 
with strongly suppressed spontaneous emission. Hence, the 
underlying mechanism is applicable to any polarizable parti- 
cle and, in particular, also to molecules where no closed exci- 
tation cycle is necessarily available. As momentum diffusion 
(heating) originates mainly from field fluctuations due to the 
randomness of photon loss via the cavity mirrors, the attain- 
able atomic temperature is limited by the cavity linewidth k 
1 3, 6]. In high-finesse cavities, one even gets sub-Doppler 
cooling (Tfl, but at the expense of a slow cooling rate. 

Among several suggestions to circumvent this slowdown, 
like e.g. using collectively enhanced light scattering 1 8], it has 
been recently proposed to add a gain medium into the cavity 
to achieve lower temperatures and fast cooling simultaneously 
The central idea here amounts to obtain an effectively 
small cavity linewidth via intracavity amplification in a not 
so good resonator In a first approach to this idea, Vuletic 
predicted improved cooling by a strong reduction of the ef- 
fective cavity linewidth through addition of intracavity gain 
j9J. In essence, the partial compensation of the cavity round 
trip loss by gain gives rise to an enhancement of the cavity-to- 
free-space scattering ratio substantially above unity, even for 
rather lossy resonators. As, in such systems, the cooling force 
is proportional to this scattering ratio, he predicted both fast 
cooling and low final atomic temperatures. Of course, gain 



in quantum mechanics is inevitably connected with fluctua- 
tions, so that the validity of the connection of temperature and 
linewidth is not so obvious any more in this case. 

In the present work, we follow this general idea but use a 
more refined approach based on microscopic modeling of gain 
by an ensemble of inverted atoms. The effects of quantum 
fluctuations and saturation are thus automatically modeled in a 
self-consistent way and the model is a direct generalization of 
our recent investigations of a single-atom-single-mode laser 
including light forces 1 10, 11]. In this configuration, we could 
show that by a suitable choice of parameters (i.e. blue detun- 
ing of the field mode with respect to the atomic transition), 
single-atom lasing can be combined with trapping and cooling 
of the atom at a cavity field mode antinode. In the stationary 
limit, the average kinetic energy of the particle can be lower 
than the optical potential depth of the laser field that is created 
by the atom itself. Extending this laser model to many atoms, 
we directly get a microscopic model for a hot cavity by simply 
fixing the positions of some of the atoms. These represent our 
gain medium as they are partially inverted through external 
driving. Of course, beyond merely providing for gain, they 
play a much more complex role in the combined atom-field 
dynamics. 

From a different point of view, we are extending the stan- 
dard Haken-Lamb laser model fl^ fisil by including light 
forces on the active atoms and study lasing as well as their 
trapping and cooling properties. As compared to the single- 
atom laser model, the trapping field is now generated from 
all the active atoms simultaneously, so that each one has to 
be less inverted to achieve a desired optical potential (photon 
number). This should reduce spontaneous emission and mo- 
mentum diffusion of the atoms and we can expect better trap- 
ping and lower temperatures. In a first step to simply mimic an 
active medium, we will keep all atoms but one fixed at antin- 
odes and study the dynamics of the remaining atom. In the 
parameter regimes where we find trapping and cooling of the 
test atom, its parameters can then be prescribed to the other 
atoms to find a self-consistent final solution for all atoms. 

Our work here is organized as follows: in SeclUwe present 
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the general model and define the relevant quantities and pa- 
rameters. In Sec.|III|we numerically study the semiclassical 
limit of the equations where we replace the field and atomic 
operators by their expectation values and factorize expecta- 
tion values of operator products. A laser model which in- 
cludes second order quantum correlations is presented and 
analytically solved in Sec. IIV Al The results are then used 
as a hot-cavity description for the motion of an extra atom in 
Sects. HVBI - IIV El Again we can find analytic expressions 
for the relevant physical quantities allowing a thorough dis- 
cussion of the properties of the present system. In particular, 
we focus on the dipole force and the optical potential for a 
fixed atom in Sec. IIV Cl In Sec. lIVE)! we calculate the fric- 
tion and diffusion coefficients which give rise to an estimate of 
the atomic localization and equilibrium temperature. Finally, 
the scaling properties of these quantities with the number of 
atoms are discussed in Sec. IIV El 

II. MODEL 



The steady state of this system so far is obviously a vacuum 
cavity field with all atoms in the ground state. Here, we are 
interested in an active system and, since gain requires atomic 
population inversion, we need to add pumping of the atoms. 
In practice, such pumping is achieved via an auxiliary atomic 
level actively coupled to the lower level which decays pre- 
dominantly into the atomic upper state within a short time. 
As introduced by Haken already some decades ago 1 13], this 
process can be modeled in a quantum mechanically consistent 
way by a so-called inverted heat bath approach which properly 
accounts for the fluctuations that are accompanied by gain, 
without the need to enlarge the atomic Hilbert space. Basi- 
cally, it amounts to introducing a spontaneous absorption rate 
2v, analogous to spontaneous emission but in the opposite di- 
rection. Mathematically, we simply have to add the following 
Liouvilhan term 

CuQ =1^^ (2(T+„igl(T_„i - {o--„iCr+,„, q}) (4) 
m 



In this section, we briefly review key elements of the 
Haken-Lamb laser model OQl, which ignores atomic mo- 
tion, to get a first insight into the dynamics on a simple and 
intuitive level. The laser active gas consists of M identical 
two-level atoms (represented by their Pauli operators (T-m, 
tJ+m, and az,m) which interact with a single mode supported 
by an optical high-finesse resonator Each atomic dipole is 
coupled to the cavity mode (represented by the bosonic field 
operator a) with coupling constant g times the mode func- 
tion f{x) evaluated at the corresponding position of the atom 
x„i. In the simplest (one-dimensional) case, the mode func- 
tion is given by f{x) — cos{kx) with mode wavenumber k. 
The cavity resonance frequency ojc — kc is generally detuned 
from the atomic transition frequency cja by a small amount 
A = — (^a, on the order of a few atomic linewidths at most. 
In a frame rotating at the cavity frequency, the atom-field 
Hamilton operator including interaction then reads within the 
rotating-wave approximation 

^int = — A a+m<^-m 
m 

-'ig^f{x„i){a'^(J-m-cr+nia)- (1) 

m 

In addition to this coherent dynamics, the atoms as well as the 
resonator mode are also weakly coupled to the environment. 
Using a standard master-equation approach fTT] for sponta- 
neous emission of the atoms into the vacuum (with rate 27 for 
each atom) and photon loss via the mirrors (with rate 2k), we 
get 

g= -i [Hint, g] + CjQ + Cf,g , (2) 
where the damping operators are given by 

^'(Q = (2o'-m6'cr+m - {cr+mCT-jn, q}) , (3a) 

ni 

CkQ = K (2aQa^ — {a^a, p}) . (3b) 



to our master equation (ME). Putting all terms together, the 
resulting ME is equivalent to the Heisenberg-Langevin equa- 
tions (HLE) [with Grn = gf{xm)], 

a = -Ka + ^ GmCT'-rn + , (5a) 

m 

&-rn = (iA - 7 - v)a^„i + GrnO-z,ma + S,a.rn , (5b) 

+2{v - 7) + C^,™ . (5c) 

Fluctuations are properly included here in the noise terms 
which follow from Eqs. Q and ©. Notice, that at this point 
the atomic center-of-mass (COM) positions are simply con- 
tained as c-number parameters il5ll . As these equations are 
coupled and nonlinear, a direct analytical solution neverthe- 
less seems out of reach and a vast amount of theoretical work 
on laser theory has been devoted to their solution in the past 
decades Iil2i] . Here, we will only resort to the most common 
proven approximations in the following. 



III. SEMICLASSICAL LASER EQUATIONS WITH LIGHT 
FORCES 

Starting from the corresponding quantum HLE for the 
polarization and field dynamics, we obtain approximate c- 
number equations by simply taking expectation values of 
these equations and factorizing expectation values of opera- 
tor products into products of single-operator expectation val- 
ues. Within this approximation, we do not need to further 
specify the noise terms as their expectation values are zero by 
construction. This approach is generally known as semiclas- 
sical laser theory. The replacement ((a), (a-m), {o'z,m)) ~^ 



3 



(a, 



i) then yields 



(iA - 7 - iy)sm + GmZma , 



In standard approaches, Eqs. (|6jl contain the positions of the 
atoms as fixed parameters via Gm — gf{xm)- As a new el- 
ement, we now introduce the atomic COM motion into the 
system so that the positions Xm [t] get dynamical and are gov- 
erned by light forces. Within the point particle model which is 
based upon a classical treatment of the atomic COM motion 
the corresponding Newtonian equations of motion 

are 

im=Pm/"T.at, (7a) 
Pm = -Fm , (7b) 

with the force acting on the mth atom given by 

F„ = -2 Im{a*s™} (V™G™) (8) 

and mat denoting the atomic mass. Note that this description 
of the external dynamics via point particles is only valid as 
long as the respective atomic kinetic energies are well above 
the recoil limit E^cc ~ ^^fc^/(2mat), which is well fulfilled 
down to atomic temperatures in the order of only a few micro 
Kelvin. 

Although the particles are now moving, the system can still 
posses a quasistationary state. Here the field as well as the 
atomic momentum and position distributions which might be 
characterized by an effective kinetic temperature are approx- 
imately constant in time. As in other single-mode-two-level 
laser models, we find a steady-state energy balance relation 
for the total atomic ground and excited state populations [i.e. 
-Pg = Em(l-^m)/2andPo = Em(l+^™)/2' respectively] 
and the mean photon number N = |a|^. 



(6a) 
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(Note that Pg + Pc = M.) This relation describes the balance 
of energy fed into the system by the pumping mechanism and 
the photon loss via atomic spontaneous emission and cavity 
field decay. Relation (|9jl is quite universal and valid beyond 
the semiclassical treatment so that we will recover it in the 
quantum model in Sects. iTVAl and lrVCl Note that it does not 
include the kinetic and potential energy fed into the atomic 
motion in single emission or absorption events (recoil energy) 
since it is only a negligible fraction of the photon energy or 
internal atomic excitation energy and, thus, only very weakly 
modifies the conditions for the lasing threshold. Thus the sys- 
tem should exhibit coherent oscillation if 



7 + 1^ 
J/ — 7 



(10) 



where we have introduced the rate of emission into the res 
onator mode per atom as 

(7 + 1^)5' 




FIG. 1: Surface plots of the common population inversion z (left) 
and the photon number A*' — \o?\ (right) as a function of the atom 
number M and the pumping rate u in the stationary state for the 
situation where all atoms are fixed at field antinodes. The parameters 

are (7, A, <?) = (!, 20, 5)«;. 



Of course, the total rate wtot = ""^Sm Pi^m) depends on 
the atomic positions and hence the spacial distribution of the 
atoms has influence on whether threshold is reached or not. 
In the limiting case where all atoms are located close to field 
nodes, the total rate wtot gets very small, which will eventu- 
ally prevent the system from reaching threshold. In this case, 
the single-atom excitation is determined by the rates of pump- 
ing and spontaneous emission, i//(7 + i/), and the cavity mode 
remains unpopulated. Hence we also get not forces on the 
atoms, and the optical potential vanishes. On the other hand, 
when the system operates above threshold, the steady-state of 
the population inversion for each atom is fixed by 



(12) 



w 



(7 + vY + A2 



(11) 



and the mean photon number can be determined from the re- 
lation (|9}. 

This behavior is depicted in Fig.n(a) where we plot the 
population inversion z as a function of the atom number M 
and the pumping strength v in the case of lowest possible 
threshold with all atoms fixed at mode antinodes (zm ~ z). 
With growing v, the population inversion reaches its maxi- 
mum value at threshold and eventually drops back to zero, 
when the growing stimulated emission (gain) decreases the 
upper state population. The cavity photon number, on the 
other hand, increases linearly with both v and M as shown 
in Fig.[T](b). 

Let us now investigate the effect of light forces on the dy- 
namics. In this context it is interesting to note that for inverted 
atoms, the effect of intensity gradients is reversed as compared 
to thermal atoms. Hence, as discussed in Refs. 1 10, 11], an in- 
verted atom is a high-field seeker for a blue detuned field. As 
a consequence, we have to work with a blue detuned cavity, 
i.e. A > 0, if we want our gain atoms to be drawn towards 
field antinodes where the gain is large. Luckily, this coincides 
with the parameter regime where the motion can be cooled via 
gain as some kinetic energy is needed to optimally fulfill the 
energy resonance condition for stimulated emission fl^. 

To demonstrate this we plot in Fig.|2](a) the trajectories ob- 
tained from a numerical simulation of Eqs. Q and Q for an 
ensemble of ten atoms where the atoms start at random posi- 
tions with random initial velocities. The atoms are decelerated 
and eventually get trapped in the vicinity of field antinodes 
where they oscillate in single wells. Inspection of Fig.|2](d) 
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FIG. 2: Sample trajectories obtained by numerically integrating 
Eqs. ^ and for 10 atoms. The parameters are the same as in 
Fig.Qwith f = 20k. 

(a) Atomic positions in units of the wavelength A. 

(b) Cavity photon number A'^. 

(c) Mean atomic distance to the nearest antinodes in units of A. 

(d) Average of the atomic velocities in units of the Doppler velocity. 
The dotted lines in (b) and (d) correspond to the results for four 
atoms. 



which shows the corresponding mean velocity truly reveals 
that the atomic motion is strongly damped at the beginning. 
After the atoms get trapped they are still cooled further, al- 
though on a much slower timescale. For comparison, we have 
added in Fig. |5](d) the curve for lowest possible atom num- 
ber to pass threshold (dotted line) which is four for the chosen 
parameters. 

Hence we see that, in principle, gain and cooling can si- 
multaneously occur and that cooling can be significantly en- 
hanced by the presence of other atoms in the same mode. For 
favorable parameters, all the atoms have accumulated at the 
antinodes of the mode after some time and the system shows 
steady lasing at the maximum intensity possible. The classical 
model, however, leaves out the fluctuations of the light forces 
and hence cannot be used to predict the final temperature of 
the atomic gas and the long term stability of this lasing pro- 
cess. A consistent description then requires to include higher 
order quantum correlations of the operators and forces as pre- 
sented in the next section. 



IV. QUANTUM DESCRIPTION OF COOLING IN A HOT 
CAVITY 



A. Active cavity model 

In this section, we include higher order quantum correc- 
tions to the laser model presented in the previous section. Let 



us start with all atoms located at the field antinodes as found 
in the long time limit above and define the collective second- 
order operators 

n =^ ^ cr+m(7-m , (13a) 

m 

S ^ {a^cr~rn + (y+mo) , (13b) 
A - ^ (aV_,„ - a+mo) , (13c) 



such that we can rewrite the interaction Hamiltonian {!) in the 
form 



H 



int 



-MAn - MgA . 



(14) 



Again, we can derive HLE for these operators compatible 
to the ME (|3 and 0, 



n ' 



n = - 2-/Ti - + 2iy (l - IT) + T 
E = - rS - AA + iMg [E, A] + T^, 
X = - r A + AE + Tx . 



(15a) 
(15b) 
(15c) 
(15d) 



It turns out that the collective operators defined in Eqs. ( fT3t . 
together with the photon number operator $ ~ a^a, almost 
form a linear, closed set of equations. The only terms leading 
out of the set are atomic cross correlations of the form 



(16) 



as well as the commutator in Eq. ( I15ck It is generally well es- 
tablished that the atom-atom correlations turn out to be small 
at least for sufficient distant atoms (d ^ A) which see un- 
corrected vacuum fluctuations 1 17]. In Eqs. (I15> we have in- 
troduced second-order noise operators that originate from the 
coupling of the system to the environment. Later we will need 
expectation values of their autocorrelation functions which are 
directly related to the damping (diffusion) terms in the ME. 

Looking at the equations for the expectation values of the 
photon number operator N — ($) and the mean upper state 
population P — (H), we recover the steady-state energy bal- 
ance equation introduced in the previous section by setting the 
time derivatives to zero and adding Eqs. ( I15a> and ( I15bk This 
relation now reads 



Mv{l-P) = kN + MjP 



(17) 



and allows to determine the steady-state photon number. 
Equations ( I15t are exact but still nonlinear as they incorpo- 
rate an operator product in Eq. ( I15c> . Due to this nonlin- 
earity, the corresponding equations for the expectation values 
are not closed and have no explicit analytic solution. Anal- 
ogous to the semiclassical case, in principle, an additional 
equation for ([E,A]) is needed which in turn will involve 
higher order moments resulting in an infinite hierarchy of 
equations. At this point we will again use a factorization ap- 
proximation for higher order products as introduced already 
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FIG. 3: (a) Steady-state photon number per atom n vs pumping rate 
V for different atom numbers M, as indicated. Tlie parameters are 
(7, A, 3) = (10,40,4)^.. 

(b) Corresponding effectively reduced linewidth HeS ■ 



in earlier laser models 0, ITsil . It relies on the replacement 

(2TI - 1)$ (2n - 1}$ such that 



iM [S,A] « 2zM* + 2n. 



(18) 



Here, zm is a scalar parameter denoting the average single- 
atom population inversion. It has to be determined self- 
consistently during the calculations of the steady state. This 
assumption will drop part of the population-intensity correla- 
tions by applying approximation ( I18> while it keeps the quan- 
tum correlations between the atomic polarization and field 
which are pronounced in lasers. 

In a next step, we make use of the small linewidth of the 
cavity with respect to the atomic decay rates, which allows us 
to adiabatically eliminate the operators 11, S and, A such that 
we remain with a single differential equation for the photon 
number operator. 



-2k(1-C)$ + 2M?7 + T<, 



(19) 



From this equation we see that the atoms introduce a gain term 
for the cavity field with an effective strength given by 



77 



^ -\- V + w 



(20) 



In this way, they compensate the round-trip loss of the res- 
onator light and induce laser oscillation. The gain parameter 
defined by 



(7 + v)i 



k(7 - 



(21) 



is proportional to the population difference zm and is posi- 
tive for positive zm (population inversion). In this case, the 
medium effectively reduces the effective resonator damping 
according to 



Kcff = k(1 - C) ■ 



(22) 



Since the cavity-induced cooling force in a passive cavity is 
inversely proportional to the loss rate, we already can expect 
the cooling efficiency to be enhanced for ^ > 0. On the other 



hand, since the atoms couple to the vacuum modes, they in- 
troduce additional noise into the system. The coiTesponding 
term reads 



T— — 



w 



"f + f 
9 



(23) 



As already mentioned, we have to calculate the value of 
the single-atom population inversion zu in a consistent way. 
Therefore, we take the expectation values of Eqs. il5\ . solve 
them in the steady state, and impose the condition 



(24) 



This yields 



ZM 



-(k(7 + ly + w) + Mw{i' — 7) 



2Mw(7 + ly) 
— |4ku'(7 + iy)M{j — v + w) 

+ [k(7 + iy + w)+ Mw{v - 7)]^}^^^) 



(25) 



Note that from Eq. \25\ it is straight forward to show that < 
1. 



B. Atomic motion in a hot cavity 

Let us now investigate the Ught forces inside the active res- 
onator by adding a single additional active atom at position 
X with the position dependent coupling G = gcos{kx). In 
some previous work, we have investigated such a configura- 
tion without the other atoms present, where the atom itself 
has to provide for the entire gain. Even then, lasing, cool- 
ingj and trapping can coexist for a blue detuned light field 
1 10] under suitable conditions. We found that atomic equilib- 
rium temperatures below the Doppler limit are possible, but 
this poses very stringent conditions onto the atomic pumping 
mechanisms. Characteristically, the final atomic equilibrium 
temperature decreases with higher intracavity photon numbers 
and the particle's localization is enhanced. 

Here we consider already Af active atoms present in the 
resonator which contribute the major part of the cavity photon 
number. Hence, we can expect to get a deeper optical po- 
tential and faster mechanical cooling rates. Of course, as in- 
evitable for any quantum consistent amplification mechanism, 
the atoms representing the active medium will introduce addi- 
tional noise in the field. This results in enlarged fluctuations of 
the system-operator expectation values, which in turn will im- 
ply also stronger momentum diffusion. The combined effect 
of these two contributions will be worked out in the following. 

As the total population inversion needed for crossing the 
laser threshold is now shared among many atoms, we get less 
stringent pumping requirements and a strongly reduced de- 
pendence of the cavity field on the individual atomic positions 
as compared to the single-atom (passive-cavity) case. In par- 
ticular, there will be still photons and thus an optical potential 
present, even when the extra atom is located very close to a 
field node. 
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Generally, we expect that larger photon numbers lead to 
lower temperatures since cooling relies on energy dissipation 
by photons that leave the resonator via its mirrors. In the ideal 
regime, each photon carries away some atomic kinetic energy. 
As a second advantage of more atoms, we expect that heat- 
ing due to spontaneous emission is diminished as the atoms 
share the required excitation. Finally, since the equilibrium 
temperature for cavity cooling scales with the linewidth k | 6|, 
an effectively gain-reduced linewidth lets us further expect en- 
hanced cooling. 

In what follows we consider a single atom moving under 
the influence of the field created by stimulated emission of the 
atom itself and the active medium described in Sec. lIV Al We 
therefore have to add the contribution of the moving atom to 
the Hamiltonian Q. Writing 11 = cr+(7_, E = aV- + cr+a, 
and A = —i{a^a^ — cr+a), with G{x) — gf{x), we have 



n = Hint - An - G{x)A . 



(26) 



Using again the factorization approximation (211 — 1)$ w 
Z$, the corresponding linearized HLE for the field and the 
extra atom read 

$ = -2Kcff $ + G(.x)S + 2AIf] + + , (27a) 

fl = -27n-G(a;)E + 2i^(l-n) + Tn, (27b) 

t = -rs - AA + 2G{x)Z^ + 2G{x)U + Ts , (27c) 

A = -FA + AS + Ta . (27d) 

These equations, containing quantum correlations up to sec- 
ond order in the field amplitude and atomic polarizations, will 
be the basis of aU our further considerations. 



C. Photon number and forces in the steady state 

As already mentioned above, as a crucial step to obtain 
an explicit approximate solution, we replace the popula- 
tion inversion operators az,m and = cr+a^ — a^a+ by 
their expectation values in the nonlinear operator products of 
Eqs. (I15c> and Mlc\ involving the photon number operator $. 
In practice we proceed as follows. We first determine zm+i 
from Eq. ( 125 > by assuming that all M + 1 atoms are fixed at 
antinodes. This population inversion zm+i then fixes the ef- 
fective linewidth n^s and we can get the population inversion 
of the moving atom by solving for (11) — [Z + l)/2 in the 
steady state. After some algebra we get 

1 



Z = 2W^(^ + ^) (^'^g(^ + u + W) + W{u--i + 2Mrj) 
- {4KeffM^(7 + v)i'j~v + W) 

+ Ksh + iy + W) + W{v - 7 + 2MT])f)'''^) . 

(28) 

Note that this approximation holds well as long as the extra 
atom is near a field antinode, where the spatial variation of 
the coupling strength is small. If the motion of the atom is 
sufficiently slow, the dynamics of the internal variables can 



adiabatically follow any atomic displacement and we can ap- 
proximate the state of the system by the steady state for an 
atom at fixed position. This allows to calculate the average 
force acting on the atom (dipole force) from the expectation 
value of the force operator, given by the spatial derivative of 
the interaction term in the Hamiltonian iflal . 



-\IU = (VG)A. 



(29) 



Since any photon in the cavity originates from stimulated 
emission, the stationary photon number N should directly be 
related to the atomic excitation. Indeed, by including the con- 
tribution of the moving atom, we find a relation analogous to 
Eq. ([T7|l, 



Mv{\ - P) + v[\ -P) = kN + M-fP + jP . 



(30) 



Notice that this equation reduces to Eq. (I17> . when the cou- 
pling of the atom to the lasing mode vanishes, i.e. at field 
nodes. Figure |3] (a) shows the stationary photon number per 
atom n — N / {M + 1) as a function of the pumping rate v for 
(7, A, g) = (10, 40, 4)k when the atom is fixed at an antin- 
ode. We plot here only the range 1/ > 7 since lower pumping 
strengths would prevent population inversion. Each curve ex- 
hibits a threshold value for v which evidently approaches 7 
with growing M. As one would further expect for our las- 
ing system, the light intensity increases linearly with v above 
threshold. Yet, the fact that n does not remain constant for 
increasing M but increases as well, demonstrates the coop- 
erative emission properties of the system. It should therefore 
be possible to achieve higher gain parameters by incorporat- 
ing many atoms for stimulated emission. This behavior is il- 
lustrated in Fig.|3](b). It depicts the corresponding effective 
resonator linewidth Kcff that is directly related to the gain pa- 
rameter C via Eq. (I22> . For > 7, the effective resonator 
linewidth is reduced linearly with v until the system passes 
threshold. There, the population inversion and hence KeS have 
already reached their minimum values and only slightly de- 
crease. As this happens close to threshold, involving many 
atoms implies enhanced gain for fixed parameters. We can 
therefore hope for improved cavity cooling since the final 
atomic temperature scales with the resonator linewidth there. 

In Fig.0](a) and (b) we show the dependence of, respec- 
tively, n and Z on the atomic position within half a wave- 
length for different atom numbers M. The change in the cor- 
responding light potential [c.f. Fig.0](d)] results in a net force 
acting on the atom which can be calculated from the steady- 
state solution of Eqs. ( I27> . It is given by 



F : 



G 



(VG) 



(31) 



where we have defined the determinant of the linear system. 



DcS = Kcff(7 + 'y + W) - + v)WZ . 



(32) 



From Eq. (13 H we see that the dipole force is an odd function 
of the atom-field detuning A. In analogy to the single-atom 
case 1 10], the atom is pushed towards regions with large inten- 
sity (high field seeker) for a blue detuned laser field. A plot of 
F for A > can be found in Fig.|3(c). 



(a) x/x (c) 1/4 




(b) x/A (d) x/A 



FIG. 4: Dependence of (a) the photon number per atom n, (b) the 
atomic population Z, (c) the force F, and (d) the hght potential V 
on the position x of a single moving atom for varying atom number 
M that constitute the gain medium. The parameters are the same as 
in Fig.|2| 

When the atom is confined to the vicinity of field antinodes, 
the spatial variations in the quantities W and Z drop. In this 
case, the force mediated by the active medium with M « M+ 
1 1 atoms, in units of the single-atom force Fq, is given by 



Fo Koff 1 - C 

This result indicates a gain-induced enhancement of the lo- 
cal forces close to field antinodes which, in general, are small 
there due to low field gradients. The increase of the dipole 
force and hence the light potential leads to strong localiza- 
tion of the atom. It should be emphasized that localiza- 
tion effects become particularly important for the velocity- 
dependent (friction) forces within the scope of the cooling 
performance. However, let us refer at this point to the dis- 
cussion at the end of Sec. lIV Dl 

D. Atomic temperature and localization 

Let us now drop the assumption of slow atoms and investi- 
gate atomic motion in more detail. In particular, we calculate 
the particle's equilibrium temperature T associated with its 
COM motion. Comparing fceT to the optical potential depth, 
we then get information about localization effects. Using the 
well known Einstein relation 1 15|, the temperature of the atom 
can be estimated by the ratio of the spatially averaged momen- 
tum diffusion term T> and the linear friction coefficient (3, 

knT^V/^. (34) 

Following the procedure of Ref. | iT], we can obtain both of 
these quantities from the force operator ( I29> . 




10 20 30 10 20 30 

M M 



FIG. 5: (a) Position averaged friction coefficient /3 in arbitrary units 
vs A/ for (7, V, g) = (10, 20, 4)k and different values of A, as indi- 
cated. 

(b) Average of the total diffusion coefficient Vtot for the same pa- 
rameters. 

(c) Approximate atomic equilibrium temperature in units of the 
Doppler temperature. 

(d) Ratio of the kinetic energy and the light potential as an indicator 
for localization of the particle. 

We will not present the whole calculation here apart from 
mentioning the basic ideas. As long as the atom moves much 
less than a wavelength within the characteristic time scale of 
the internal dynamics, the friction can be approximated by the 
force term linear in the particle's velocity v. A simultaneous 
expansion of both Eqs. i27\ and the system operators in terms 
of V results in a set of dynamical equations, where the first 
order solution yields the friction coefficient 

/3 = (VG)(Ai). (35) 

The full expression of (A^) is listed in App.lAl 

For the diffusion coefficient we have to evaluate the two- 
time correlation function of the force operator LI 9 1 since mo- 
mentum diffusion arises from fluctuations of this operator 
Here we will use the approach of Ref. |20] as it allows us 
to directly relate momentum diffusion to the noise caused by 
spontaneous emission and cavity decay. We therefore read off 
V from 

{F{t)F{t - r)) - {F{t)){F{t -t))= 2V5{t) . (36) 

A very brief derivation can be found in App.|B] We also have 
to include heating from the random recoil of spontaneously 
emitted photons, which results in an additional diffusion term 
that can be found in Ref. flsll . 

In Fig.|5l(a) and (b) we depict the scaling of, respectively, 
the position averaged friction and total diffusion coefficients 
with the amount of atoms in the active medium, for differ- 
ent values of the detuning A. (Note that we plot continuous 
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curves instead of discrete values here, and that M = corre- 
sponds to the one-atom case.) We chose v — 2^^ in order to 
ensure population inversion such that (3 is negative. Clearly, 
we see that both (3 and I?tot strongly depend on the number 
of atoms, and their absolute values grow with M. Since each 
additional atom increases the dipole force F and thus the light 
potential V [see Fig.|3(c) and (d)], one could expect this ef- 
fect to arise from enlarged light intensities when many atoms 
constitute the medium. Indeed, due to the fact that friction re- 
lies on the time-delayed force caused by an atomic motional 
displacement, we should conclude that (3 shows similar scal- 
ing properties with M as the dipole force F. On the other 
hand, the fluctuating part of the force operator quantifies mo- 
mentum diffusion, not its steady-state value. As we saw from 
Eq. \\9\ . the active medium introduces an extra noise term 
which should lead to additional momentum diffusion. More- 
over, we infer from Fig.|5](a) and (b) that larger values of the 
atom-field detuning A, certainly decreasing the photon num- 
ber, result in larger friction as well as momentum diffusion. 
However, we will later see how we can compare situations 
with different atom numbers in a more favorable way where 
the photon number, the atomic population inversion, and the 
gain parameter are independent of AI. Let us finally mention 
here that the cooling rate is inversely proportional to (3 such 
that large values of (3 will lead to fast cooling. 

The resulting equilibrium temperature T associated with 
the atomic COM motion is shown in Fig. |5] (c) in units of 
the Doppler temperature Td = hj/kB- Analyzing this plot, 
we can identify two regimes according to the scaled pumping 
rate y ~ i^/A. As long as there are only a few atoms in the 
cavity, the atomic temperature decreases with increasing y. In 
fact, this is in strong analogy to the single-atom case where 
we showed that temperatures even below the limit of passive 
cavity cooling are possible for sufficiently large y fllll . Here, 
on the other hand, the addition of even a single atom can lead 
to a significant drop in the temperature, which already proves 
the importance of cooperative aspects. Depending on the ac- 
tual value of the detuning A, we eventually find, with growing 
M, a transition point to a second regime where the cooling 
performance can be enhanced for lower pumping rates y that 
we achieved by raising A. The cooling enhancement only 
happens for sufficiently large atom numbers and is directly re- 
lated with the fact that strong pumping immediately saturates 
the gain here. Indeed, one has to ensure not to saturate the gain 
as the cooperative effects of emission will then break and, as 
a consequence, the cooling efficiency is reduced. Each curve 
shows a minimum which tends towards larger M when we 
increase the detuning A. For fixed parameters, the minimum 
corresponds to a somehow optimum atom number M which 
keeps the balance between the required population inversion 
on the one hand and saturation on the other hand. In this sense, 
lowering the scaled pumping rate y means that we can involve 
more atoms into the cooling process without saturating the 
gain at once. Note that for the chosen parameters here, which 
do not require to be in the strong atom-field coupling regime, 
one can achieve sub-Doppler cooling for M > 20. 

In view of a setup where one has to assure that the atoms 
remain in the interaction region close to field antinodes, we 



will now focus on localization effects. Therefore we plot in 
Fig. m (d) the ratio of the atomic kinetic energy E and the 
light potential V. Large values of this ratio correspond to a 
more or less even distribution of the particle's position while 
values below one indicate strong localization. One can again 
identify the afore mentioned regimes of low and high atom 
numbers M, respectively. Qualitatively, the curves for E/V 
are similar to the ones for the temperature and we recognize 
that cooling is directly accompanied by atomic trapping for 
sufficiently high M. In particular, the ratio E/V drops well 
below unity when the atomic temperature gets down into the 
range of the Doppler limit. 

As already mentioned before, the cooling properties may 
significantly depend on the initial atomic temperature and lo- 
calization. For rather cold and well localized atoms, in par- 
ticular, the position dependence of the cooling forces must be 
taken into account. Typically, the cooling time scale strongly 
increases as soon as an initially rather hot ensemble of cavity- 
cooled atoms gets trapped. (Notice, that we have already 
seen a similar effect in Sec.|lIl]although no information about 
the equilibrium temperature and cooling time could be given 
there.) Moreover, it has been shown that in the limit of very 
strong confinement, the friction forces can even vanish such 
that the cooling process completely stops I^J). On the other 
hand, strong cooling forces may arise in the high-intensity 
regime resulting in very strong atomic localization 
where the additional momentum diffusion leads to a mere 
slight increase of the atomic equilibrium temperature. 

In order to get an estimate of the influence of localization 
effects on the cooling performance in a hot cavity, we will de- 
rive the gain-induced enhancement of the friction for a rather 
well localized atom. For 7\f « Af + 1 1, it is given by 

where (3o indicates the single-atom friction coefficient. In the 
previous section, we have calculated the same expression for 
die stationary force which turned out to scale inversely with 
the effective linewidth Kcff- Since (3/Po i^cS^ cool- 
ing forces are even more enhanced than the trapping potential 
such that an atom trapped very close to field antinodes still 
experience relatively large friction as soon as 1 — ^ ^ 1. 
We therefore deduce that, besides providing strong particle 
localization, the stimulated-emission mediated cooling forces 
can significantly accelerate the cooling process as compared 
to conventional schemes using passive resonators. 



E. Scaling properties witli the number of atoms 

In the good-cavity limit (k ^ v), the gain-induced cool- 
ing of a single particle can lead to atomic temperatures well 
below the limit of passive cavity cooling for sufficiently large 
amount of pumping ratio y and when one can neglect spon- 
taneous emission 1 11]. Nonetheless, the required large exci- 
tation rate i/ certainly poses a serious experimental challenge 
there. As we already have seen from Fig.|5j the use of many 
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beneath the single-atom hmit y = 1/2. The corresponding ra- 
tio E/V is depicted in Fig.|5](b). Again we find simultaneous 
trapping and cooling. This result suggests the possibility to 
use a rather bad cavity and many weakly pumped active atoms 
as gain medium to achieve stimulated emission induced cool- 
ing. Fortunately, these conditions simultaneously ease the re- 
quirements on cavity technology and on the necessary amount 
of pumping per atom y and give a rather promising perspec- 
tive of possible implementations of the scheme. 



FIG. 6: (a) Equilibrium temperature T in units of the minimum tem- 
perature achievable by passive cavity cooling. T is plotted as a func- 
tion of the atom number AI where k/{M + 1) was kept constant for 
different pumping rates y, as indicated, 
(b) Corresponding ratio E/V. 

atoms for stimulated photon emission into the resonator can 
enhance trapping and also lower the particle's kinetic energy, 
even for small values of y. In the following, we will now 
analyze the scaling of the cooling properties with respect to 
the atom number M to see whether this method should work 
for larger ensembles under less restrictive experimental con- 
ditions. 

For this purpose, we introduce a rescaling of the system 
parameters as follows. As we have seen above, the compari- 
son of the motional properties for different M is actually not 
trivial. Particularly, not only the photon number TV but also 
the population inversion zm as well as the gain parameter 
play a crucial role there as they all depend on M in a highly 
nonlinear way. A good choice will certainly be a rescaling 
which keeps these three quantities constant at all. A closer 
look at Eqs. ( I25> , ( I2H . and il9i reveals that one can realize 
this all at once by keeping the quotient k/ (A/ + 1) constant. 
Consequently, the operating point of the lasing system can be 
set independently from the atom number without any further 
complications. We will accomplish this here by relating the 
total amount of emission into the resonator to the photon loss 
rate according to 

{M+l)w = aK (38) 

via the parameter a. The value of a and hence the operating 
point can then be set by a proper choice of the coupling con- 
stant g. 

Besides, this rescaling accounts for comparable conditions 
also in the energy balance. Having fixed a certain operating 
point, the established light intensity no longer depends on the 
number of atoms involved. Regardless, however, the rate of 
photon leakage through the cavity mirrors, given by 2kN, lin- 
early grows now with the atom number Since each photon 
carries away energy and entropy from the system, we expect 
not only lower temperatures but also faster cooling. 

In Fig.|6l(a) we show the scaling of the atomic temperature 
with M for different values of y where we assumed 7 = 0. 
Clearly, T drops well below the limit of passive cavity cool- 
ing which is given by k-gTcav — fiK. When many atoms are 
involved, this happens even for values of the pumping ratio y 



V. CONCLUSIONS 

We have shown that intracavity gain significantly modifies 
light forces in cavities. In particular, the stimulated-emission 
dominated interaction of individual particles with a single 
light mode opens new perpectives on the dynamics of lasing 
and light-induced manipulation of the atomic motion. To ex- 
hibit the underlying physics in detail, we developed a quan- 
tum consistent theory for intracavity gain, based on a micro- 
scopic description of an ensemble of individual inverted atoms 
strongly coupled to a single resonator mode. Including light 
forces into well-established laser equations, we found a self- 
consistent analytical model for the combined dynamics of the 
atomic motion in a hot cavity. As an extra bonus, by account- 
ing for spontaneous emission and cavity decay, fluctuations 
and saturation are automatically included in our model by con- 
struction. 

As a central result, we can confirm previous claims that 
laser cooling can be strongly improved in such a setup with 
less stringend requirements on cavity technology. This result 
turns out to be closely related to a strong reduction of the 
effective resonator linewidth as the amplifying atoms partly 
balance the photon loss via stimulated emission. In this way, 
photon scattering into the cavity, which accounts for the cool- 
ing forces, can be drastically enhanced. While the naive ex- 
pectation of a much reduced temperature due to the narrow 
linewidth is not fully met by the results, we still get much 
faster cooling and sub-Doppler temperatures. Our findings 
show the potential usefulness of stimulated-emission induced 
cooling in optical resonators and may open the way to novel 
experiments. 

More importantly however, we recover that optical lasing 
can be concurrent with strong atomic localization and sub- 
Doppler cooling. This directly points to the possibility of 
combined photon lasing and atom lasing in the continuous- 
wave (CW) regime. In a conceivable joint setup, a coherent 
atomic beam could be out-coupled from an ultra cold ensem- 
ble of atoms in an optical lattice. The atoms are cooled to the 
motional ground state without particle loss by Raman gain in- 
volving stimulated emission into a resonator mode. The pho- 
ton leakage through the cavity mirrors provides the photon 
laser channel there. In order to achieve combined CW lasing 
operation, the cooling region is fed by atomic tunneling from 
neighbouring sites. Ground-state atoms which are prepared in 
a non coupling Hyperfine state by a proper microwave beam 
leave the trap due to gravity. 
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APPENDIX A: FRICTION FORCE 



Evaluating the correlation function ( IB1> . one finds for the con- 
tribution of the active medium to the diffusion coefficient 



(7 + I/)2A2G2Z2(VG) 



Mw 



r(7 + t/ + w)2 
X {(7 + vf{^) + [k(7 + vf + (7 - u)wT] (TT) 
+ (7 + i^)r/.g [k(7 + V + w) + w{j — i^)] 

X - /r2(S) +,y[wT+{^ + l^f] } (B3) 



We have 



(Ai) = - 



Aw^{yz) 



[vW + M77(7 + + W)] 



+ (7 + z.)2(r-2«eff)^]} 
AW^(VM^) 



[k^sv + M?7(7 + v)Z] 



X (2Kcff(7 + v)T[k^s{i + V-W) 
+ (7 + u)WZ] - G'jKK^ff - (7 + y? 
x{T-K,s)Z + W[n,s-{l + v)Zf]). (Al) 

APPENDIX B: DIFFUSION COEFFICIENT 

Writing F = (F) + T, Eq. (|36j becomes 

(T(t)T(i - r)} - 22?5(t) , (Bl) 

where the fluctuating part of the force operator reads for quasi- 
stationary conditions and for D — (F^ + A2)Deff 

'^-^{(7 + ^)AGZ(T$ + T^) 

+ Keff AGTn + Kcff (7 + y) (ATs + FTa) 

+ G2[At,ff-(7 + z.)Z]TA}. (B2) 



by collecting all terms involving the operator T-^. Here we 
used the steady-state solution of the collective operators 



(n) 



^ + v + w 
2u)i^+(7 + ^)zM+i($) 
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^ + V + w 



The remaining portion finally reads 



(B4a) 
(B4b) 
(B4c) 



V^,,,^^^^^{n{^ + vf/\^G'Z^[vW 

+ Mr;(7 + v^W)\+ KeffA^G^ 

X [k^sW + 27(Kcff - WZ) + Mriiu - -f)WZ] 

+ 2VI^A2 [Keffi/ + Mr7(7 + v)Z] 

X [Koff(7 -y) + «(7 + y)Z] 

X {2«:eff (7 + i^) + ['^off - (7 + y)Z] GVr} 

+ M-q [(7 + i/)(7 + 1/ + T4^) - kWZ]] 

X {(F^ + A2)«;2^(7 + + G^ [/teff - (7 + ^^)^]' 

+ 2«:eff (7 + '^)rG2['*cff - (7 + v)Z])) . (B5) 
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